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Abstract: Following the analysis of [1], we calculate the spectrum of fluctuations of a probe 
DA:-brane in the background of Dp-branes, for k = p,p -\- 2,p + A and p < 5. The result 
corresponds to the mesonic spectrum of a (p + l)-dimensional super- Yang-Mills (SYM) theory 
coupled to 'dynamical quarks', i.e., fields in the fundamental representation - the latter are 
confined to a defect for k = p and p + 2. We find a universal behaviour where the spectrum 
is discrete and the mesons are deeply bound. The mass gap and spectrum are set by the 
scale M ~ ruq/ geffimq), where ruq is the mass of the fundamental fields and geff{mq) is the 
effective coupling evaluated at the quark mass, i.e., giffijrnq) = gyj^^Nrnq ^. We consider 
the evolution of the meson spectra into the far infrared of three-dimensional SYM, where the 
gravity dual lifts to M-theory. We also argue that the mass scale appearing in the meson 
spectra is dictated by holography. 
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1. Introduction 

The AdS / CFT correspondence and its extensions have provided a powerful framework for the 
study of strongly coupled gauge theories in various dimensions. The original correspondence 
[2, 3], that certain conformal field theories are equivalent to string theory on AdS backgrounds, 
was extended in [4] to a more general gauge/gravity duality in various dimensions. Itzhaki 
et al [4] studied the general case of a stack of A^" coincident Dj»-branes in the limit in which 
brane modes decouple from the bulk. They argued that the super- Yang-Mills U{N) gauge 
theory on the {p + l)-dimensional worldvolume of the Dp-branes is dual to the closed string 
theory on the 'near-horizon' background induced by the brancs. 

These dualities follow from a straightforward extension of the usual decoupling limit [5] 
now applied to Dp-branes [4] . For general p (7^ 3) , the gauge theory is distinguished from the 
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conformal case by the fact that the Yang-Mills coupling gvM is dimensionful. Hence there is 
a power-law running of the effective coupling with the energy scale U : 

9lff = 9lMNUP-^ . (1.1) 

The duality relates the energy scale and the radial coordinate r transverse to the Dp-brane 
worldvolume in the usual way, U = r /a' . In the dual background, the absence of conformal 
invariance in the general case is manifest in the radial variation of both the string coupling 
(or dilaton) and the spacetime curvature — see below for details. As the supergravity back- 
ground is only trustworthy for weak string coupling and small curvatures, it provides a dual 
description of the theory which is reliable for an intermediate regime of energies. In this 
regime, the dual gauge theory is always strongly coupled. 

In a complementary direction, Karch and Katz [6] demonstrated that probe D7-branes can 
be used to introduce fundamental matter fields into the standard AdS / CFT correspondence 
from D3-branes. Inserting n D7-branes into the AdS5 x background corresponds to coupling 
n flavours of 'dynamical quarks' {i.e., n hypermultiplets in the fundamental representation) 
to the original four-dimensional SYM theory. Adding these extra branes/ficlds also reduces 
the number of conserved supercharges from sixteen to eight. The hypermultiplets arise from 
the lightest modes of strings stretching between the D7- and D3-branes and their mass is 
niq = Ljl-Kol where L is the coordinate distance between the two sets of branes (and as 
usual, Xjl-Kol is the string tension). The resulting gauge theory containing quarks has a rich 
spectrum of quark-antiquark bound states, which henceforth we refer to as 'mesons' [1]. In 
the decoupling limit, the duals are open strings attached to the D7-branes and the calculation 
of the meson spectrum in the field theory becomes an exercise in studying the fiuctuation of 
probe branes. These ideas have been further developed in a number of directions towards the 
goal of constructing gauge/gravity duals for a QCD-like theories [7, 8, 9] and in particular, 
the meson spectrum has been studied in a number of different contexts [1, 10]. 

In this paper, we use the gauge/gravity duality to explore the meson spectra of such 
SYM theories containing fundamental fields in different numbers of spacetime dimensions. 
In particular, following [1], we calculate the spectrum of fluctuations of a probe D/c-brane 
super symmetrically embedded in the background of N Dp-branes, for k = p,p + 2,p + 4 
and p < 5. This corresponds to the mesonic spectrum of a C/ (A^) super- Yang-Mills theory in 
p + 1 dimensions coupled to a hypermultiplet in the fundamental representation. Again, these 
dynamical quarks arise as the lightest modes of the {k,p) and {p,k) with their mass given 
by rnq = L/2Tra' . Also as before, the number of conserved supercharges is reduced to eight 
by the addition of the fundamental hypermultiplet. A consequence of the supersymmetric 
embedding of the probe brane is that the quarks are conflned to a defect of codimension two 
and one fov k = p and p + 2, respectively. So it is only for A; = p -|- 4 that the D-brane 
configuration yields fundamental fields propagating in the full p+1 dimensions of the gauge 
theory.^ 

^However, for the k — p and p + 2 cases, one might consider compactifying the Dp-brane worldvolume 
directions transverse to the Dfc-brane, as in [7]. 
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The resulting meson spectra for all of these different configurations display certain uni- 
versal characteristics, which are common to the original D3-D7 results [1]. In general, the 
spectra are discrete and the mesons are deeply bound. Up to numerical coefficients, the meson 
masses can all be written in terms of a single scale: 

M ^mq/geff{mg) , (1.2) 

where is the quark mass and geff{mq) is the effective coupling (1.1) evaluated at this 
mass scale, g1^^{mq) = ^fyjy^iVmq" . In particular then, eq. (1.2) gives the mass gap of the 
spectrum. The detailed calculation of these results is presented in section 2. 

As an example, for N background D2-branes and a D6-brane probe, we find that the 
meson masses scale as M oc rriq'^ j (^g^ j^^N^^'^ ■ It is interesting to note then that in [11], a 
similar study found M oc for the same configuration of branes. The resolution of this 
apparent discrepancy is that the latter analysis considers very low quark masses in the far 
infrared, i.e., niq <^ gyM, where the dual gravity configuration actually lifts to M-theory. To 
better understand this behaviour in different regimes of the field theory, in section 3, we study 
probe D2- and D4-branes in the D2-brane background, for which computations can easily be 
extrapolated between the type IIA and M-theory regimes. Hence, in this section, we calculate 
the meson spectra in (2 -|- l)-dimensional SYM coupled to a fundamental hypermultiplet 
confined to codimension-two and -one defects, in the far infrared regime of the theory. 

Finally, section 4 presents a discussion of our results and future directions. In particular, 
we present a simple argument that the mass scale of the meson spectra is dictated by the 
consistency of holography and wc briefly describe several tests of the latter idea. A number of 
technical results arc provided in subsequent appendices. While this paper was in preparation, 
we became aware of [12] which addresses the same basic problem as considered here, however, 
we provide some complementary analyses and a different interpretation of the results. 

2. Supergravity meson spectra in various dimensions 

In this section, we study the excitations of a DA;-brane probe (with k = p + 4,p + 2,p) 
supersymmetrically embedded in the 'near horizon' geometry induced by N coincident Dp- 
branes (with p < 5).^ These calculations yield the spectrum of mesonic states in the dual 
(p + l)-dimensional field theory with eight supercharges. We first review the background 
geometry induced by the Dp-branes and then proceed to compute the spectra of fiuctuations 
of the probe D/c-branes. 

2.1 Background geometry 

The supergravity solution corresponding to N coincident Dp-branes is, in the string frame 
(see, e.g., [14] and references therein) 

^Thoro is no dual field theory for p > 6 as no decoupling limit is possible [4, 13]. We comment on the 
interesting case of p = 5 in section 4. 
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^2^ = /;3-p)/2 Co..., = -(/-^ - 1) (2.1) 

where ds^(E'^^'P-') denotes a fiat spacetime metric with one time and p spatial directions. 
The Y^, with A = \, ...,9 — p, parametrize the (9 — p)-dimensional space transverse to the 
Dp-branes.'^ The harmonic function depends on the transverse radial coordinate r = \Y\ 

/, = ! + (?)"' M 

where the constant rp is defined in terms of the number A'^ of Dp-branes, the string coupling 
constant Qs, and the inverse string tension a': 

rl'P = g,N (47ra')(^-^')/2 r(2=£) /An . (2.3) 

Following [4], we take the decoupling limit (in which the open string modes propagating on 
the Dp-branes decouple from the bulk closed string modes) with 

9ym = (2vr)^-25. a'^^"'^/' = fixed, a' ^ 0. (2.4) 
This limit also holds r/a' constant [2] so that eq. (2.1) reduces to the 'near horizon' solution: 

d,^ = (^)""'"%»^(Ea,P)) + (i)"->/V.<if 

This supergravity background then provides a dual description of U (N) super- Yang- Mills 
theories in p+1 dimensions, the worldvolume field theory on the Dp-branes. In accord with the 
duality, both the background and the field theory have sixteen supersymmetries. The isometry 
group for the background geometry (2.5) induced by the Dp-branes is SO{l,p) x 50(9 — p) 
[4]. From the perspective of the dual gauge theory, SO{l,p) corresponds to the spacetime 
Lorentz symmetry while 50(9 — p) is the R-symmetry group. 

As mentioned in the introduction, the supergravity sohition (2.5) is a trustworthy back- 
ground provided that both the curvatures and string coupling arc small. This limits the 
supergravity description to an intermediate regime of energies in the field theory or of ra- 
dial distances in the background. In terms of the effective coupling (1.1), this restriction is 
succinctly expressed as [4] 

1 < geff < . (2.6) 

Hence the field theory is strongly coupled where the dual supergravity description is valid. For 
p < 3, the spacetime curvature is large for very large values of r, invalidating the supergravity 
description. However, in this UV regime, the coupling runs to geff <C 1 and so perturbative 



^Thcsc directions correspond to x"* with A = p + 1, p + 2, • • • , 9 in the arrays presented below — see, e.g.; 
eq. (2.7). 
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field theory is applicable. For small values of r, both the effective coupling and the dilaton are 
large. However, another dual description can be found to describe this infrared regime with 

4 

9eff ^ -/V^-p [4]. For p > 3, the coupling runs in the reverse direction. Hence perturbative 
SYM applies in the far infrared where one finds large spacetime curvatures for small r. The 
effective coupling grows in the UV and a new diial theory must be found for very large r. 
For p = 3, the theory is conformal and the supergravity background (2.5) can be used for all 
energy regimes [2]. 

2.2 Meson spectrum in Dp-D(p + 4) 

Consider the following configuration of N coincident Dp-branes (0 < p < 5) and one D(p + 4)- 
brane probe, in which the D(p+4)-brane is parallel to the Dp-branes' worldvolume directions: 

1 •••pp+lp + 2p + 3p + 4p + 5- -- 9 
Dp X X ••• X (2.7) 
D{p + 4)xx---x X X X X 

Embedding the D(p+4)-brane in the (p+5) • • • 9 plane at |y| = L, the configuration of branes 
remains supersymmetric but only half of the original sixteen supercharges are preserved. 
Correspondingly, the SO{% — p) symmetry of the background geometry (2.5) is broken to 
SO{A) ~ SU{2)r X SU{2)l acting on yiy2y3y4 for L 7^ (which is the case of interest 
here) SOiA — p) acting in the remaining transverse space around the D(p + 4)-brane. The 
R-symmetry group for fields on the D{p + 4)-brane is SU{2)ji x S'0(4 — p). The fields can be 
classified according to their transformation properties under SU{2)ji x SU{2)l x >S'0(4 — p).^ 
The induced metric on the D(p + 4)-brane probe is, from (2.5), 

7— p 7— p 

ds' = ' "^''^"^^'^'^^^[t^Tw] ' (^^^' + ^'^^3) (2.8) 

where = r"^ — L? and are spherical coordinates in the (p + 1) • • • (p + 4)-space. The 
dual gauge theory is {p + l)-dimensional U{N) SYM coupled to a hypermultiplet of matter 
fields in the fundamental representation. The quark mass is the distance L separating the 
D(p + 4)-brane from the D^?-branes multiplied by the string tension l/27rQ;': 

^ (2.9) 
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We wish to calculate the spectrum of mesons corresponding to open string excitations on 
the D(p + 4)-brane — the analysis closely follows that of [1]. These mesons are represented 
by excitations of the scalar and gauge worldvolume fields. Their dynamics is governed by 
the D(p + 4)-brane action. The relevant terms of the latter include the Dirac-Born-Infeld 



*We note, however, that the 50(4 — p) properties tend to be less useful. In particular, this group only 
really acts for p < 2. 
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where 



(DBI) action and only the term with C(p_|_i) from the Wess-Zumino action (see, e.g., [14] and 
references therein): 

S = Sdbi - j P[C(j,+i)] AFAF (2.10) 

Sdbi = -rp+4 J d^'+Ve"^ V- det(P[G]„6 + 27ra'F„b) . (2.11) 

Here, rp+4 = 27r/(7s(47r^a;')(*'+^)/^ is the D(p + 4)-brane tension and, as usual, P denotes the 
pullback of a bulk field to the probe brane's worldvolume [15]. The ten-dimensional spacetime 
metric Gab and the bulk Ramond-Ramond (RR) form C(p+i) were given in (2.5).^ 

Note that the action (2.10) represents the bosonic part of the full action invariant under 
eight supercharges and that it has SU{2)ji x SU{2)l x SO{A — p) symmetry corresponding to 
rotations in the transverse space. The Wcss-Zumino term breaks the symmetry interchanging 
the SU{2)r and SU{2)l. In the dual gauge theory, this corresponds to the asymmetry of 
SU{2)l and SU{2)r: the former commutes with the supercharges while the latter, as the 
R-symmetry group of the theory, does not. 

We begin by considering fluctuations in the position of the probe D(^?-|-4)-brane. Working 
in the static gauge, these are scalar fluctuations about the fiducial embedding: 

= L5^ + 2'Kax'^, ^ = 5,.. .,9 -p. 

We only need to retain terms in eq. (2.10) to quadratic order in these fluctuations and so the 
relevant part of the Lagrangian density for fluctuations in the scalar fields is 

(27ra')2 /rp\(7-p)/2 



(2.12) 



where gab denotes the induced metric (2.8) on the D(p + 4)-brane. Summation over the 
repeated index A is implied. The factor e~'^^J— det gab is independent of the fluctuations 
which is a reflection of the supersymmetry of the brane configuration. The latter dictates 
there be no potential for the position of the probe brane. Retaining terms only to quadratic 
order in the fluctuations, we can drop terms containing from the factor g'^ jr^~v)l'^ and 
so we can use (2.8) as the induced metric on the probe brane in the final result. 

The equations of motion for each of the 5 — p fluctuations resulting from the variation of 
the Lagrangian density (2.12) are: 



da 



(p2+L2)(7-p)/4' 



0, (2.13) 



■''The index labelling conventions here and in the following arc: The indices a, b, c, . . . denote the worldvolume 
directions of the probe brane. Greek indices n, v denote probe brane directions parallel to the background 
branes. For directions in the probe brane worldvolume orthogonal to the background branes, we use spherical 
polar coordinates with radius p and angular coordinate indices denoted by indices i, j, k. Finally, indices 
A,B,C, . . . denote background directions orthogonal to the Dp-brane worldvolume. 
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where we have taken x to be any one of the fluctuations. Also hij is the metric on the 
unit three-sphere which, along with radial coordinate p, spans the (y-^,...,y^) directions. 
Expanding the equation of motion, we obtain: 

j^^j^^,d,d^X + ysMP%x) + ^Mai X = (2.14) 

where is the Laplacian on the unit three-sphere. Next, using separation of variables, we 
write the modes as 

X = 4>{p)e"'-''y''{S') , (2.15) 
where {S^) are spherical harmonics on the satisfying 

v^^y'' = -hih + 2)y'' , (2.16) 

and transforming in the (^, ^) representation of 50(4) = SU(2)l x SU{2)r. Then, substi- 
tuting (2.15) into (2.14), and setting 

„ U2 7-P 

^^ = -^' (2-17) 
we obtain the following equation for (f){p) = 4'{g): 

dim + + ( , ttLv. - ^^^^1 <t>i^) = • (2-18) 



(1 + p2)(7-p)/2 g2 

The solution to this differential equation must be real- valued and regular [1]. The solutions 
must also be normalizable in order to be dual to a meson state in the field theory. Thus, 
real solutions were chosen that were regular at the origin. The eigenvalues M were then 
determined by requiring that the solutions were convergent at p — > oo. For p = 3, the 
solutions can be determined in terms of hypergeometric functions [1] - see Appendix A.^ 
However, for general p, there are no known analytic solutions to this equation and one must 
resort to numerics. 

Solving (2.18) to find the eigenfunctions (P{q) and the dimensionless eigenvalue M, yields 
the (p + l)-dimensional mass spectrum of mesons mass 
eigenvalues will depend on the radial quantum number n (which corresponds to the number 

of nodes in (f){g), the radial profile) and the angular momentum quantum number £3. Hence 
we denote the eigenvalues for these scalar fluctuations as M'^{n,l^). 

Using (2.4) and (2.9), the spectrum of excitations can be expressed in terms of field 
theory quantities as follows: 



= 4^ ( ) m2 . (2.19) 



5-p / 2P-27r(P+i)/2 ' 



^Analytic solutions can also be found for eq. (2.18) in the case p = 5 but one finds there are no normalizable 
solutions! 
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Further then from eq. (1.1), we see that the meson masses scale a M oc niq/ geffiniq). 

Turning to the fluctuations of the gauge fields on the probe brane, the equations of motion 
for these fields follow from (2.10) as 

da (e-V-detffedi^'^^) - p (p^ + L^^i^-P)/^ ^u, ^ j^. ^ g ^ (2.20) 

where is a antisymmetric tensor density on the three-sphere, taking values ±1. Also V[3]j 
is the covariant derivative on the of unit radius, and coordinates are indexed as described 
in footnote 5. The first term comes from the Dirac-Born-Infeld part of the action, while the 
second is from the Wess-Zumino term and is only present if h corresponds to an index. 

We proceed to solve (2.20) by expanding and Ap in (scalar) spherical harmonics 
and Ai in vector spherical harmonics. There are three classes of vector spherical harmonics. 
The first is the covariant derivative on the three-sphere of the usual scalar spherical har- 
monics, V[3]j3^^='(S'^), while the other two, labelled [Vj^^'^(-S'^), have ^3 > 1, transform in the 
(^, ^) of 50(4), and satisfy 

VLJ^f = (2.21) 



where Rj = 26j is the Ricci tensor on the three-sphere of unit radius. 

We first consider modes satisfying d^A^^ = 0. In this case, the equation of motion for A/j, 
decouples from the other equations of motion, so we have the following two types of modes: 

Typel: = C^'/'iIpK''" -V'^ ^ • C = 0, ^Ip = 0, A, = (2.22) 

and 

Type 2: ^^ = 0, A^ = Mp^''"' (S^) , A = MpV"-' ^ [3]iy'' {S^) ■ (2.23) 

Modes containing the yf^'"^ spherical harmonics will not mix with other modes because they 
are in different representations of 50(4) [1]. Hence a third independent set of modes is 

Type 3: A^ = 0, Ap = 0, yl, = 0±(p)e*" '±(53). (2.24) 

There may be other modes with d^A^ 7^ 0. Such modes with /c^ = need not be 
considered because they do not yield regular solutions. On the other hand, modes with 
d^Afji 7^ and A;^ 7^ can always be put in a gauge so that they become type 2 modes and 
hence these modes are equivalent to modes discussed above in the d^^A^ = gauge. 

We now compute the spectrum for each type of mode. For type 1 modes, the equations 
with b = Q and b = i are automatically satisfied, so we need only consider the equation with 
b = jji. This equation simplifies from (2.20) to 

r^^^ 1 1 

+ -^dp{p^dpA^,) + -^v^l = , 
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and, upon substitution of (2.22), becomes 

^d,{g%MQ)) + ^^-^^^M9) - ^-^^^Ue) = 0. (2.25) 

The {p + l)-dimensional mass spectrum of mesons is then given by (2.19) with the dimen- 
sionless eigenvalues M = A{[i{n,i^) determined by this equation. However, expanding out the 
first term in eq. (2.25), one finds that it precisely matches the previous radial equation (2.18) 
for the scalar fluctuations. Hence the spectra of these two sets of fluctuations are identical, 
i.e., Mi{n,is) = MsinJs). 

For type 2 modes with £3 = 0, the gauge fields on the three sphere vanish, i.e., = 0. 
The only nontrivial solution of cq. (2.20) (setting b = /i) is then cl)2{p) ~ p~^- However, 
requiring a regular solution at the origin forces us to take the trivial solution ^2(p) = 0. Thus 
the physical type 2 modes only exist for £3 > 1. 

Putting b = jj in (2.20) and substituting (2.23) for the modes, the type 2 radial modes 
(4 > 1) satisfy 

-^dp {p^Mp)) = 4(4 + 2)Mp)- (2.26) 

Using (2.26), the equations obtained from (2.20) with b = p and b = p are equivalent. After 
using (2.17), we obtain 

d, Q d.ie'MQ))) - 4(4 + 2)4>2iQ) + ^i+^2^i-p)/2 MQ) = o, (2.27) 

and the solutions of this equation determine the spectrum of meson masses which we denote 
as M|(n, ^3). In this case, if we expand out cq. (2.27) in terms of (/){g) = q4'2{q)j we find that 
the result again matches eq. (2.18). Hence the spectrum of the type 2 modes precisely matches 
that of the scalar fluctuations and the type 1 gauge fluctuations with M2(?7.,4) = Ms{n,£3) 
with 4 > 1. 

Finally, for type 3 modes, the equations (2.20) with h = p, and b = p are automatically 
satisfled. The equation with 6 = i, an coordinate, becomes: 

P \ Tp ^ J p^Tp 

-(I^(p2 + L2)(^-^)/^e,,,V[3],A, = 0. 

Tp ^ 

Substituting (2.24) for Aj, using the identities (2.21), and making the substitutions (2.17), 
we find the following equation for (f)f{g): 

7-ri 

n2^ 



T(7-p)(4 + l)(l + ^')^0^(^) = 
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Once again, the spectrum of meson masses is given by (2.19) with the values of M = 
Ma -l- (71,^3) determined from this equation. The spectrum of mesons can again be related 
to the other spectra with M3^±{n,is) = Ms{n,is ± 1) with £3 > 1. While in general the 
spectra can only be evaluated numerically, the previous relation is established analytically 
by mapping eq. (2.18) to cq. (2.28). One finds that the type 3,=b modes with £3 = L can be 
related to the solutions (P{q) of eq. (2.18) with £3 = L ±1 yiaJ 

^3,es=L = Q% [q'-%s=l-i] (2.28) 

There is actually a subtlety for the ^3 profiles. Typically the radial ODE's yield asymp- 
totic solutions which behave like (f){g) = big~'^^ + b2Q~°''^ in the limit ^ — 00. The mass 
eigenvalues are then determined by requiring that the coefficient bi vanishes for the term with 
ai negative. However, in certain cases with (p'^ , eq. (2.28) yields solutions where both ai are 
positive. Specifically for this mode, we have cti = 6 — p — £3 and 02 = ^3 + 1 and so for 
£3 < 6 — p, both modes converge at infinity. However, in these cases, supersymmetry 'selects' 
the physical profile as that with ^3 + 1, i.e., this choice yields a super symmetric spectrum. 

2.2.1 Analysis of the spectrum 

In the previous subsection, we computed the spectra of bosonic mesons in (p+ l)-dimensional 
U (N) SYM coupled to a fundamental liypcrmultiplet by computing the spectrum of fluctua- 
tions of scalar and gauge fields on the D{p + 4)-brane worldvolume. Classifying the massive 
meson states in representations (ji,j2, js) of SU(2)jiX SU{2)lX 30(4— p), where ji,2,3 denote 
the SU{2)ji^L and SO{4 — p) spin, the bosonic modes of the D(p -|- 4)-brane give rise to the 
following mesonic states and mass spectra:^ 

• A — p scalars in the (^, ^, 1) corresponding to fluctuations transverse to the Dp-branes, 
with mass Ms{n,£s), n > 0, £3 > 0; 

• 1 scalar in the (^,-^,0) corresponding to transverse fluctuations but parallel to the 
separation of the Dp- and D(p-|-4)-branes, with mass (72,^3), n > 0, ^3 > 0; 

• 1 vector in the (^, y^^) corresponding to type 1 gauge fields on the D(p -|- 4)-brane, 
with mass Mi{n,£3) = MsinJs), n > 0, ^3 > 0; 

• 1 scalar in the (y-, ■^)0) corresponding to type 2 gauge fields on the D(p + 4)-brane, 
with mass AhinJs) = Ms{n,£3), n > 0, ^ > 1; 

^In more detail, the relation (2.28) is established as follows: Set = L — 1 in eq. (2.18). Re-express the 
equation in terms of = Q^~^F{q), where the power was chosen to eliminate the term (L^ — 1)10^ 4>. Next 
one multiplies the result by (1 + j)^)^"?')/^ and takes a derivative with respect to q. F only appears in the 
resulting equation through the derivatives F', F" and F'" {i.e., there are no terms proportional to F) and 
so the result can be re-expressed as a second-order equation in terms of G{g) = F'{g)/g^, which one finds 
matches precisely cq. (2.28) for (f)'^ ^^^i^. 

*Keep in mind that the 5*0(4 — p) group only really acts for p <2. 
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• 1 scalar in the (^%-'-, 0) corresponding to type 3, + gauge fields on the D(p + 4)- 
brane, with mass M3^+(n, ^3) = Ms{n,i3 + 1), n > 0, ^3 > 1; 

• 1 scalar in the (^^5^, ^^^^,0) corresponding to type 3, — gauge fields on the D(p + 4)- 
brane, with mass M3^_(n,£3) = (71,^3 — 1), n > 0, £3 > 1. 

In summary, the spectra can be related to each other through 

M,(n,4) = Mi(n,4) = M2(n,4) = M3,+ (n,4 - 1) = M3,-(n,4 + 1) ■ (2.29) 

Again, for general p, we were unable to find an analytic solution for the spectrum but the 
mass eigenvalues are fixed by a simple ODE (2.18) and may be determined numerically. An 
analytic solution was found in [1] for p=3 - see Appendix A. Another special case is p = 1 
where analytic radial profiles (written in terms of Bessel functions) can be found for £3 = 0. 
As advertised in the introduction, all of the masses scale parametrically as M oc ruq/ g^.^ f{mq), 
where geff{mq) is the running effective coupling (1.1) evaluated at the quark mass scale. 

As noted earlier, the Dp-D(p + 4) brane system preserves eight supercharges which corre- 
sponds to M = 2 supersymmetry in four dimensions. The mesons, as massive representations 
of the supersymmetry algebra, should fill out long supermultiplets. The construction of these 
supermultiplets for the {p + l)-dimcnsional theories with eight supercharges is a simple ex- 
tension of that for M = 2 multiplets in the four-dimensional theory, considered in [1]. The 
multiplets are constructed by acting with the supercharges Q on a state with spin ^ under 
SU{2)r which is annihilated by the Q's. Since SU{2)l commutes with the supercharges, all 
states in a given supermultiplet will be in the same representation of SU{2)i^. Of course, as 
dictated by supersymmetry, each multiplet contains an equal number of bosonic and fcrmionic 
components, i.e., 8(^3 + 1) of each. For £3 > 2, the bosonic contributions to the generic mul- 
tiplet are: 6 — p real scalars and one vector in the ^ of SU (2)j^ and two real scalars in the 
^ lb 1 of SU{2)ii. Now the fermions can be represented as Dirac spinors with two {p = 1, 2) 
or four (p = 3, 4) components. Hence, for p = 3, 4, the generic multiplet contains a single 
Dirac fermion, in each of ^^±1 and ^^zii of SU{2)fi. For p = 1,2, there arc two Dirac fermions 
in each of the ^^^^ and . However, in the latter case, there is also a nontrivial 5*0(4 —p) 
action and these pairs of spinors combine together in the spin-half representation of this ro- 
tation group. ^ Exceptional semi-short supermultiplets appear for £3 = 0, 1. For £3 = 0, the 
spectrum contains 5 —p scalars and one vector which are singlets of SU{2)ji, one scalar in 
the 1, and one (p=3,4) or two (p=l,2) Dirac fermions in the ^. For £3 = 1, there are 6 — p 
scalars and one vector in the i, one scalar in the | and one (p=3,4) or two (p=l,2) Dirac 
fermions in each of the and 1. 

In each of the cases above, the bosonic content is, of course, identical to that found using 
supergravity. Hence supersymmetry allows us to extend the computed bosonic spectrum to 
include fermions. First we have for p = 1,2 

®These considerations may be extended to p = 0, where the dual theory is supersymmetric matrix quantum 

mechanics coupled to degrees of freedom in the fundamental representation. The latter include four fermionic 
variables, which are organized as a spinor of the 5*0(4) symmetry. 
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• 1 fermion in the (^2±i, ^) with mass M/,i(n, ^3) = Ms(n, ^3), n > 0, £3 > 0; 

• 1 fermion in the (§, ^) with mass M/,2(n,4) = Ms(n,t^ + 1), n > 0, £3 > 0. 



Keeping in mind that the spinor representations arc twice as large for p = 3, 4, we can keep 
the classification as above except for dropping as there is no 50(4 — p). 

2.2.2 Meson spectrum for p = 2 

As a specific example, let us consider the case p = 2 in which the background geometry is 
induced by N coincident D2-branes and there is one D6-brane probe embedded a distance L 
in the 789-directions: 

0123456789 
D2 X X X (2.30) 
D6 X X X X X X X 

The radial differential equation (2.18) (for scalar fluctuations of the D6-brane) is 



3 / /\f2 



m = 



(2.31) 



which can be solved using the shooting method. For either ^ ^ or ^3 ^ 00, the solution has 
the form ^(^) = Ag^^ + Bg~^^^^. Using the boundary condition (f){Q ^ 0) = g^'^, we solved 
(2.31) numerically. As the solutions to (2.31) must be regular for all values of g, we tuned 
the constant M to obtain the regular solutions, which behave as (?!) = g~^^~^ for ^ — 00. In 
this way, the three-dimensional mass spectrum of scalar mesons was found to be 



471" ,-r2 

— M 

3 -2 



m" 



9ymN 



(2.32) 



where the values of M5 (71,^3) are given in table 1. As noted earlier, the eigenvalues for the 
type 1 and 2 gauge fields are identical to those for the scalar modes and hence their spectrum 

is also given by (2.32). For the type 3 gauge field modes, we solved (2.28) using the shooting 
method and found eigenvalues identical to those in table 1 but with the £^ label shifted as 
indicated in eq. (2.29), i.e., M3,±(n,£3) = Ms{n,i3 ± 1), £3 > 1. 





n 







1 


2 


3 


4 


5 







11.34 


36.53 


75.49 


128.19 


194.65 


274.86 


is 


1 


33.39 


70.37 


121.01 


185.36 


263.41 


355.20 




2 


66.20 


114.96 


177.33 


253.36 


343.08 


446.49 




3 


109.75 


170.30 


244.43 


332.17 


433.57 


548.65 



Table 1: The mesonic spectrum in terms of M^(n,£3) corresponding to a configuration a D6-brane 
probe in the D2-brane geometry. Each mass eigenvalue is labelled by the quantum numbers £3 and n 
(where n represents the number of nodes in the (j){p) function). 
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2.3 Meson spectrum in Dp-D(p + 2) 

Consider tlie following configuration of A'^ coincident Dp-branes (1 < p < 5) and one D(p + 2)- 
brane probe, embedded a distance L from the D]?-branes: 

l---p-lpp + lp + 2p + 3p + 4- -- 9 
Dp X X ••• X X (2.33) 
D{p + 2)xx---x X X X 

As before, this orientation of branes was chosen to produce a supersymmetric system which 
is in static equilibrium. Inserting the D(p + 2)-brane breaks the SO{l,p) symmetry of the 
01 • • • p-directions to SO{l,p— 1) acting on 01 • • • (p— 1). Similarly the SO{9—p) symmetry of 
the space transverse to the Dp-branes to SO{3) = SU{2)ji acting on yiy2y3 ^^^^^ f^j, L ^ 0, 
SO {5 — p) acting on • ■ ■ Y^~'p. The R-symmetry group for fields on the D(p + 2)-brane is 
SU{2)ji. The induced metric on the D(p + 2)-brane is 

7— p 7—p 

2 



ds' = [^^^j ' + (tJ^) ' (dp' + P'dnl) (2.34) 

where, as usual, we have introduced spherical polar coordinates (with radial coordinate p) in 
the probe brane worldvolume directions orthogonal to the background Dp-brancs. 

The dual gauge theory is {p + l)-dimensional, but the fundamental hypermultiplet has 
been introduced on a p-dimensional surface. Hence the matter fields are localized on a 
codimension-one defect [16], e.g., = 0. The theory again has eight conserved supercharges 
and the quark mass is still given by (2.9). Together the rotation symmetries above, i.e., 
SU{2)ii X SO{5 — p), form the R-symmetry group of the gauge theory. 

We follow the same procedure described for the Dp-D{p + 4) system to compute the 
spectrum of mesons corresponding to fluctuations of the D(p-|-2)-brane. In this case, however, 
the D{p + 2)-brane can fluctuate in the X^'-direction, parallel to the background branes, 
as well as in the ■ ■ ■ y^~*'-directions, transverse to the background Dp-branes. We take 
X"^ {A = 4, ...,9 — p) and as the fluctuations in directions transverse and parallel to the 
background Dp-branes, respectively: 

= 5^ L + 2TTa'x^, A = 4,...,9-p (2.35) 
XP = 27raV. (2.36) 

As with the Dp-D(p + 4)-brane configuration, the relevant action for the D(p + 2)-brane 
fields is the DBI action (2.11) combined with the Wess-Zumino involving C(p_)_i): 

S = Sdbi - 27ra'rp+2 j P[C(p+i)] A F . (2.37) 

In this case, the Wess-Zumino term produces a coupling between the gauge fields and the 
scalar il). 
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Considering first the fluctuations directions orthogonal to the background branes, 

the relevant quadratic Lagrangian density is 

(27ra'^2 



(2.38) 



where summation over A = A, ...,9 — p is implied. Here we have dropped the 'constant' term 
proportional to the determinant of the induced metric gab and the dilaton, since it is again 
independent of the fluctuations. Now in an expansion to quadratic order the entire pre-factor 
in eq. (2.38) can be treated as though it is independent of x^- particular, we can take gab 
to be given by (2.34). 

As each of the fluctuations appears on equal footing in the Lagrangian (2.38), the 
superscript is dropped in the following. Now we expand the modes as 



X = '^±(p)e*"y^(S'^) 



(2.39) 



where e"^^'^ is a plane wave in the 0, 1, ...,p — 1 space and y^^{S'^) are spherical harmonics 
on an S'^ of unit radius (£2 = 0, 1,2, ...). With g and M as defined in (2.17), the equation 
determining the radial profile is 

M2 



Ui2 + I] 



(1 + ^2)(7-p)/2 



Mq) = 0. 



(2.40) 



As discussed in section 2.2, the solutions corresponding to physical mesons are real, regular 
at the origin and convergent asymptotically. Analytic solutions to this equation foi p = 3 
(in terms of hypergeometric functions) are given in appendix A.^^ For p = 1,2,4, we solved 
(2.40) numerically (using the shooting method) to determine the mass eigenvalues. As in 
the previous section then, for 1 < < 4, the spectrum of mesons (here confined to a {p- 
l)-dimensonal defect) is given by (2.19) with the dimensionless eigenvalues M = M_l (71,^2) 
computed from (2.40), where n is the number of nodes of the (t>±{p) function. 

The relevant Lagrangian density for fluctuations of the gauge fields and tp, the fluctuation 
of the D{p + 2)-brane along the XP direction, follows from (2.37) as 



£ = -(2W)Vp+2 e-'^^- del gab U " ^ci^dd^J + -F^dF' 



7-p 
2 



1 



Tied 



7-p 



{Fo^dpip + F^pdei/j + Fped^tp) 



, (2.41) 



where p, 9, ip are spherical polar coordinates in the y^y^y^-space. To quadratic order, the 
equations of motion are 



da 



Vdet h ( — 



7-p 

2 



g^'dbi^ 



{7-p) 



pr 



5— p 



J-P 



■e*^V[2].A, 







(2.42) 



^"Eq. (2.40) also becomes a hypergeometric equation for p = 5 but there are no solutions satisfying all of 

the necessary criteria. 
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da VditTlF-^^) - (7 - p)p!:^e''^5iV = (2.43) 

where we are using the index notation given in footnote 5. Also hij and e*-^ are the metric 
and antisymmetric tensor density on the unit two-sphere, respectively. Note that e^"^ = 1. 
The second term in each equation results from the Wess-Zumino term and in (2.43) this term 
is only present if b is an index on the S^. 

We can expand the scalar and gauge fields in terms of spherical harmonics on the two- 
sphere component of the D{p + 2)-brane: A^,Ap,tp in terms of scalar spherical harmonics 
and Ai, in terms of vector spherical harmonics. Note that the equations of motion (2.42) and 
(2.43) imply that the scalar field ip couples only to the Ai gauge field modes. Thus, working 
in the gauge d^Afj, = 0, we can define one type of mode not coupled to as 

Typel: A^ = ^^Mpy''-^y''iS2) , ^■k = 0, A, = , A = . (2.44) 

For '0 = 0) there is a second type of mode given by 

Type 2: A^ = , = 02(p)e*"y^(52) , A, = ^2(p)e*"V[2],y^(52) . (2.45) 

Finally, there is a third type of gauge field mode which is coupled to -i/j for ^2 > 1: 

Type 3: ^^, = 0, Ap = , A^ = Mp)e"'-''Vd^eijV\2^y^'{S^) , (2.46) 

where the factor \/det h makes up for the density weight of the two-dimensional e-symbol. 

We now proceed to compute the spectra for each type of mode. For type 1 modes, the 
equation of motion is (2.18) with h = ji. Substituting (2.44) and making the redefinitions 
(2.17), we obtain 

dl4>i{Q) + -dMQ) + - M^^) Mq) = 0. (2.47) 

^ \(i + ) 2 s J 

Note that this result is identical to cq. (2.40), the differential equation for the scalar fluctu- 
ations transverse to the background branes. Thus, the mass spectrum M = Mi (71,^2) here 
will be identical to that for the transverse scalars. 

For the type 2 gauge fields and "0 = Oj (2.42) is identically satisfied. Note that for £2 = 0, 
Ai = 0. Then (2.43) with h = jj, yields (f)2 ~ p~^- As the latter is not regular at the origin, 
the only solution for £2 = is trivial {i.e., (f)2 = 0) and we need only consider £2 > 1 for these 
modes. 

For £2 > 1, (2.43) with b = 11 gives dp{p^(p2) = hih + 1)02- With this, the equations 
obtained from (2.43) with b = p and b = i are equivalent and give, with the definitions as in 
eq. (2.17), 

dlMQ) + -dsMQ) + ^Ir^ + Mq) =0, (2.48) 

& + Q^) S J 
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which defines the mass spectrum which we denote as Mll^^-^s), ^2 > 1- Putting 02 = ^(^2, 
this equation becomes the ODE (2.40) and thus the spectrum is again identical to that for 
the transverse scalars and type 1 gauge fields: M2(n, ^2) = M±{n,i2) = M\{n,l2)- 

As noted above, the type 3 gauge field modes are coupled to the scalar field tp which 
represents fluctuations of the probe branc along the direction, parallel to the background 
Dp-branes. The mode ^2 = is an exception since the gauge field vanishes: = 0. Then, 
with if) = e^^'^ (j)\\{p) and the redefinitions (2.17), eq. (2.42) yields 

dlhis) + {- + ^^-^) d,4>\\{e) + -^^hi^) = 0- (2-49) 

Solving this equation and imposing regularity requirements yields the spectrum of mesons 

M||(n,0) =M^(n,l). 

For ^2 ^ 1; we proceed via separation of variables, expanding the scalar field as ■)/' = 
e^^'^ (l)\\{p)y^'^ {S"^) and using (2.46) for the gauge field, so that (2.42) becomes 

(7)'"'^V||(P) + {p^'-^dMp)) - ^^^^ [Mp) + (7-P)^^3(p)) =0, 

(2.50) 

while (2.43) with b = i gives 

{^y~'M^Mp) + {r'-'dMp)) - ^^^^ (Up) + (7 -P)^0||(P)) = . 

(2.51) 

In both of these equations, r = \J + L'^. We diagonalize this system of equations by defining 
two new radial functions 0+ (p) , 4>- (p) as 

4>+ip) = ^Mp) + pMp) ' ^-(P) = (^2 + 1)Mp) - PhiP) ■ (2-52) 

With these new functions and also the definitions (2.17), (2.50) and (2.51), the decoupled 
equations become 

dlM9) + ^d,MQ) + (-—r^ - ^ - (7 -p)^) Me) = 

\{i+Q^)^ / (2.53) 

dlMe) + ^^d,Ms) + [-^^^^ - + (^-P)w) Me) = o. 

As usual, by solving these equations and imposing regularity requirements, the eigenvalues 
M are found for 1 < jo < 4 (while for p = 5 there are no normalizable modes). We 
denote the eigenvalues for 4>+{p) and (f'-ip) respectively as M-|_(n, £2) and M-{n, £2)- The 
spectrum of mesons corresponding to fluctuations of the D(p + 2)-brane in the p-direction 

^^Here again, one finds a subtlety for the 4>- profiles. In the limit q —> 00, the solutions behave like 
(f>{Q) = biQ~°^^ + b2Q~"^ where ai = 6 — p — £2 and 02 = £2- Hence both a, are positive for < 6 — p and 
both modes converge at infinity. Hence, in these cases, wo still determine the physical masses by demanding 
that 61 vanish, which yields a supersymmetric spectrum. 
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and fluctuations of the gauge fields is then given by (2.19) with these values of and these 
can be related to the spectrum of scalar fluctuations transverse to the background branes 
via M±{n,i2) = M±{n,£2 i 1), ^2 > 1- Again, the spectra are only evaluated numerically in 
general, however, this matching of the spectra is established analytically by mapping eq. (2.40) 
to eq. (2.53). One finds that the <p± modes are related to the <p± modes as follows: 

4>~M=L = g'-d, [^^-^0±,^,=L-i] (2.54) 

In computing the spectra of fluctuations of scalar and gauge fields on the D(j) + 2)-brane 
worldvolume, we have found the spectra of mesons living on a codimension-one defect in 
the {p + l)-dimensional super- Yang-Mills theory. Prom the supergravity computations, there 
are 6 — p scalar mesons corresponding to transverse fluctuations of the Dp-branes with mass 
M_L, one vector meson corresponding to type 1 gauge fields with mass Mi, one scalar meson 
corresponding to type 2 gauge fields with mass M2, and two scalars, corresponding to the (p± 
modes with masses Mj-. The spectra for these different modes are simply related with 

Ml{n,£2) = Miin,^) = Mi{n,i2) = M+(n,£2 - 1) = M-(n,£2 + 1) ■ (2.55) 

As with the Dp-D{p + 4) brane configuration, we were unable to find an analytic solution 
for the spectrum for general p. However, the mass eigenvalues are fixed by the ODE (2.40) 
and can be easily computed using numerical techniques. The notable exception is the case 
p = 3 for which solutions were found in terms of hyper geometric functions - see Appendix 
A. Hence while the precise numbers change, the masses scale as M oc mq/geff(mq), which is 
identical to the scaling found in the Dp-D{p + 4) brane system. 

As discussed in section 2.2.1, the mesons should fill massive supermultiplets. Counting 
the bosonic contributions found using supergravity, we see that each multiplet contains eight 
bosonic degrees of freedom, as expected from the dual gauge theory. 

2.4 Meson spectrum in Dp-Dp 

Consider the following configuration of N coincident Dp-branes (2 < p < 5) and one Dp-brane 
probe, in which the probe brane is a distance L from the background branes: 

1 •••p-2p-lpp+lp + 2p + 3- -- 9 
background x x • • • x x x (2.56) 
probe X X • • • X x x 

The orientation of D-branes was again chosen to preserve supersymmetry and the branes 
are in static equilibrium. Embedding the probe brane in this way reduces the number of 
supercharges from sixteen to eight. Correspondingly, the isometry groups of the background 
geometry (2.5) are broken. The 5*0(1, p) symmetry of the 01 • • • p-directions has been reduced 
to 5*0(1, p — 2) X 5*0(2). Further the 50(9 — p) symmetry corresponding to rotations in the 
. . . y9-P-directions is broken to 50(2)' acting in the (p + l,p + 2)-plane and for L 7^ 0, 
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5*0(6 — p) acting in the remaining transverse directions orthogonal to the separation of the 
branes. The induced metric gab on the probe brane is 



7-p 7—p 
2 / \ 2 



[^^^^j ' '^^'(^^'''"'^^ + (t^^) ' idp" + p"d^")' (2-57) 

where, as usual, we are using polar coordinates p, 9 in the probe brane worldvolume directions 
transverse to the background Dp-branes. 

The dual description has {p + l)-dimensional super- Yang- Mills coupled to a fundamental 
hypermultiplet confined to a (p — l)-dimensional surface. That is, the matter fields live on 
a codimension-two defect [16]. The quark mass (2.9) again corresponds to the mass of a 
fundamental string stretching between the background branes and the probe brane. The 
theory has eight conserved supercharges and the R-symmetry group has two components: 
the SO{Q — p) rotations and the diagonal rotations in 50(2) x 50(2)'. 

For this brane configuration, the fluctuations of the probe brane fall into two classes: 
those orthogonal and parallel to the background branes. We write the fluctuations around 
the fiducial embedding as 

= 5^L + 2TTa'x^, ^ = 3,...,9-p; (2.58) 
= 2Tra'^p^, B = p-l,p. (2.59) 

The action defining the dynamics of the Dp-brane probe is given by the DBI action (2.11) 
plus the Wess-Zumino term with (7(^+1): 

5 = Sdbi - Tp I P[0(p+i)] . (2.60) 
For the 'orthogonal' fluctuations, the quadratic Lagrangian density is 

C = -^^r.e-'^ V'^dit^ (^) f%x^dax^ (2.61) 

where summation over A is implied and = + p^. As before, the quadratic Lagrangian 
depends only on the derivatives of x^- Taking x to be any one of the x"^, we expand the 
modes as 

X = <^±(p)e*'■V^^ (2.62) 

where e**^'^ is a plane wave in the (0, 1, ...,p — 2)-space and £ = 0, ±1, ±2, .... The equation of 
motion reduces to determining the radial profile (t)±{Q)'- 

dlMe) + -d,M6) + [jYT^yfZpfJ^ ^) M9) = (2.63) 

where we use the definitions (2.17). The solutions with i = are not normalizable and 
so we restrict to £ nonzero. Since the equation above is symmetric in i, modes with the 
same absolute value of i have the same mass. For p = 3 analytic solutions in terms of 
hypergeomctric functions are given in appendix A while for p = 2,4 solutions were found 
numerically.^^ Thus, for 2 < p < 4, th e spectrum of mesons living on the (p — l)-dimensional 
^■^Here and for the remaining modes, one finds there are no normalizable modes for p = 5. 
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defect is given by (2.19) with M = M_L(n, {£ / 0) computed from eq. (2.63). 

We now turn to the 'parallel' fluctuations il^P~^,ijjP. In this case, the Wess-Zumino term 
introduces a coupling between these two fields. The quadratic Lagrangian density following 
from eq. (2.60) is 



(27ra' 



7-f 
f \ 2 



7-p 



(2.64) 

where = + again. The tpP~^,ij;P mixing can be diagonalized by working with the field 
^ = tpP^^ + itpP and its complex conjugate = tpP~^ — iipP, in terms of which eq. (2.64) 
becomes 



7-f 



7-p 



(2.65) 



and the resulting equation of motion is 



On 



r-p 

2 



g'^'db^ 



5— p 



(7-p) r_ 
2^ 'rl-^ 



= . 



(2.66) 



Proceeding via separation of variables, we expand the field as ^ = (p\\{p) e e and eq. (2.66) 
reduces to 



.(1 + ^2)V 



f (7 - p)£ 

1 + ^2 



</>ll(e)=0, (2.67) 



where we have substituted in the definitions (2.17). As before, we found analytic solutions 
for p = 3 (see appendix A) and numerical solutions for p = 2, 4. 

As seen in previous sections, there is a mapping from eq. (2.63) to eq. (2.67). One finds 
that the modes can be related to the 'orthogonal' modes via^^ 



(2.68) 



U=-L = g'-^Op [^'-^<^x,|^|=L-i] , L>2. 



Hence the spectra can be matched analytically, even if they are only evaluated numerically. 
We denote the spectra M^(n, \£\) for £ > and £ < —2, respectively. Thus it follows from 
eq. (2.68) that M+{n,£) = M±{n, \£\ + 1) and M-{n,£) = M±{n, \£\ - 1). 

The spectrum for is obtained by noting that the radial ODE for this field is identical 
to eq. (2.67) with the replacement i —£. Hence for we have: M<,^-|-(n, £) = M^(n, —£) = 
Mj_(n, 1^1 - 1) with £>2- and M^-{n,£) = M+{n, -£) = M±{n, \£\ + 1) with £<0. 



^^Note that is a subtlety here in that there is no mapping for L = —1! Hence supersymmetry seems to 
dictate that the £ = —1 solutions of eq. (2.67) are unphysical. 
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Finally, we turn to computing the spectrum for gauge field fluctuations. The linearized 
equation of motion comes entirely from the DBI part of the action (2.60) and is 

5a(e-V-det5cdi^'^^) = 5„(pF«'') = . (2.69) 

Following the earlier analysis, we find that there are two types of modes: 

Typel: = C^</)i(/9)e*V^^ k ■ ( = 0, Ap = 0, Ae = (2.70) 

Type 2: A^ = 0, Ap = Mpy'^'^'e'^^ ^ = ^2(p)e*V^^ (2.71) 

where£ = 0,±l,±2,.... 

For type 1 modes, eq. (2.69) with 6 = /x gives 

1 / /^X 

OIMq) + -dsMQ) + -JEF - ^ Mq) = 0, (2.72) 

^ \(1 + ) 2 / 

where we have used the redefinitions (2.17). This equation is identical to that for the trans- 
verse scalars (2.63) and thus the type 1 gauge fields will have the same spectrum as those 
scalar modes: Mi{n, = M±(n, \£\), £ 0. 

For modes of type 2, eq. (2.69) with b = gives 

pdp{p(P2) = -iidp^2. (2.73) 

There is no regular solution for ^ = and we only consider ^ 7^ in the following analysis. In 
this case, using (2.73), the equations resulting from (2.69) with b = p and b = 6 are equivalent 
and, using eq. (2.17), give 

3 / _ 1 \ 

OIMq) + -dMs) + -JEF - Me) = 0. (2.74) 

As usual, for 2 < p < 4, the spectrum of these fluctuations follows from solving for M = 
M2{n,£) and using (2.19). These can be matched with the spectrum of 'orthogonal' fluc- 
tuations with M2{n, \£\) = M±{n, \£\) for £ / 0. The latter follows since upon substituting 
<f>2 = i^/e, eq. (2.74) reduces to eq. (2.63). 

Thus, the spectrum of mesons in {p + l)-dimensional (2 < p < 4) SYM coupled to a 
fundamental hypermultiplet confined to a co dimension- two defect is discrete with mass gap 
M ~ m,g/geff. There are (7 — p) scalar mesons corresponding to the 'orthogonal' fluctuations 
of the probe, two corresponding to 'parallel' fluctuations, one vector corresponding to type 1 
gauge fields, and one scalar corresponding to type 2 gauge fields. The spectra of all of these 
are related via: 

M^{n, \£\) = M^{n, \£\) = A'hin, \£\) = M_{n, \£\ + 1) = M+{n,£- 1) 

= M,,+ (n,£ + l) = M,,_(n,|£| -1). (2.75) 



-20- 



As usual, the mesons organize themselves into massive supermultiplets and the masses of the 
fermionic mesons match those of the bosonic modes determined here. 

Again we have found that the meson masses scale as M oc mq/geff{mq), which is identical 
to the scaling found in the previous two cases. Examining the meson spectra of the various 
brane configurations in more detail, we note that they arc related as the radial equations for 
each different probe can be mapped into one another. Eq. (2.63) for the Dp-Dp brane system 
reduces to eq. (2.18) for the Dp-D(p-|-4) system upon substituting: (p±{Q; n,£) = g(j){Q; n,£s = 
\£\ — 1). Hence the mass levels for these two systems are identical, although the degeneracies 
will differ in the two cases. Similarly, eq. (2.40) for the Dp-D(p + 2) system maps to eq. (2.18) 
for the Dp-D(p + 4) system upon substituting: (t)±{g;n,£2) = y/Q(t^ {^Q't^^lz ~ ^2 — Of 
course, this is only a formal identification because the spectra in each case are only evaluated 
for integer values of the angular quantum numbers £2,3- Hence the physical spectra differ for 
these two systems. 



3. Beyond lOD supergravity: an example 

It is well understood that, with the exception of p = 3, the supergravity regime discussed 
above is only applicable for a well-defined intermediate regime of energy scales [4]. From 
the results of section 2, the meson spectra clearly focus on a very precise energy scale in the 
field theory, namely, the quark mass. Thus, to study meson spectra in gauge theories in the 
infrared (IR) or ultraviolet (UV) regimes with this holographic framework, one must venture 
beyond ten-dimensional supergravity. 

In this section we consider on the meson spectra beyond this supergravity regime for 
the specific example of the D2-brane background, in which we introduce D4- and D2-brane 
probes. The corresponding field theory is three-dimensional U{N) super- Yang-Mills (SYM) 
coupled to a fundamental hypermultiplet on a codimension-one and -two defects (for the 
D4- and D2-probcs, respectively) with eight supercharges. For this theory, the Yang-Mills 
coupling constant and the dimensionless effective coupling constant are given by [4] 

9'yM = ^, 9lff{mq) = '^. (3.1) 



The meson spectrum for the D2-brane theory has different descriptions depending on 
the energy scale niq [4], as shown in figure 1. Perturbative SYM is valid for g^ff <^ 1 
which corresponds to the UV regime or very large values of mq. The type IIA supergravity 
description takes over in the regime 1 < geff < N'^/^, where both the curvature and string 
coupling of the gravity background are small. Once g^ff > N"^/^ , the string coupling becomes 
large and an eleven-dimensional supergravity description is required. Thus, in the far infrared, 
the gravity theory is strongly coupled type IIA supergravity which lifts to M-theory as an 
M2-brane background. 

The choice of background here was motivated by the desire to compare our results to those 
of ref. [11]. There an extension of the gauge/gravity duality beyond the probe approximation 
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M2-branc near- core Uplifted D2-brane IIA D2-brane Perturb ative SYM 

\ \ \ * 

Sym 



Figure 1: Descriptions of the D2-brane theory for different energy regimes in the field theory, set by 
the quark mass niq. 

with a large number of fundamental fields was discussed for a system of D2- and D6-branes. 
Their results seem to indicate that the meson mass is directly proportional to the quark mass, 

3 /2 

M oc niq. In contrast, our results in section 2.2.2 gave M oc niq . These disparate scalings 
are reconciled by noting that the results of [11] actually only apply in the far infrared regime, 
where the dual description is in terms of the M2-brane throat rather than the D2 background. 
We study the theory with D4- and D2-brane probes instead of D6-branes, because it is easy 
to follow the physics of the former probes from the ten-dimensional regime to the eleven- 
dimensional phase where the latter lift to M5- and M2-brane probes, respectively. Our results 
in the previous sections indicate that the mass gap in meson spectra of these different theories 
should all scale in the same way. 

3.1 Strong coupling: M2-branes with M5- and M2-probes 

With p = 2, eq. (2.1) shows that the dilaton diverges as r — > and so in the far infrared the 
system becomes strongly coupled. As discussed above, the type IIA theory is lifted to M- 
theory in this regime. In this case, the array of D2- and D4-branes in (2.33) or of D2-branes 
in (2.56) would be lifted to a system of A'^ coincident M2-branes with an M5- or M2-brane 
probe, as indicated below: 

0123456789 11 
M2 X X X 

M5 X X X X X X 
(probe) M2 x x x 

The type IIA background solution (2.1) is readily lifted to eleven dimensions as (see, for 
example, [14] and references therein) 

= ^-2/3(_^^2 ^ ^^2 ^ ^^2^ ^ f/^{dr'^ + r^dnl + dxl^) (3.3) 
^012 = (3.4) 

where r is a radial coordinate in the 3456789-space and xn is the compact coordinate in the 
eleventh dimension with xu ~ xu + 27ri?ii. Also / = /(r, xn) is, in general, a harmonic 
function of all of the transverse coordinates. Lifting the type IIA solution (2.5) with p = 2 
does not yield exactly the M2-brane solution but rather the solution for a set of M2-branes 
smeared over the circle direction. However, we need to consider the solution for A'^ coincident 
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M2-branes localized at a point on the xn-circle. (Of course, even more complicated solutions 
can also be constructed.) In this case, the harmonic function is given by [5] 



n=—oo '■ ^ ' ' 

The summation can be carried out to yield a closed-form result — see appendix B — but 
this expression is not very illuminating for the following. 

For the fluctuation analysis, we focus on the limits r 2> Ru and r, xn <C Ru, in which the 
harmonic function (3.5) simplifies considerably. In the large r limit, r S> -Ru, the dependence 
on xii drops out leaving: 

Rii r° 

In the following, we refer to this as the "uplifted D2-brane" solution because this corresponds 
to precisely the eleven-dimensional lift of eq. (2.5) with p = 2. Near the core of the M2-branes 
{r,xii <C i?ii), the harmonic function reduces to 

/(-'-n) = p^^. (3.7) 

The metric (3.3) with this harmonic function gives the near-horizon geometry for N coincident 
M2-branes, i.e., AdS4^ x 5^ with units of flux. We will refer to this as the "near core" 
solution. 

In the configuration (3.2), we embed the probe branes at a distance L from the background 
M2-branes in, say, the X^-direction. As the system is supersymmetric, there is no potential 
for either of the probe branes and they are in static equilibrium. Working in the static gauge, 
we consider small fluctuations in the 6789-directions transverse to the background M2-branes 

= 6^L + 2Tra'x^, ^ = 6,7,8,9 (3.8) 

which are overall transverse directions for either probe brane. For the D2-brane probe, the 
rcii-fluctuations would also fall into this class. For these scalar fields, the relevant probe- 
brane action is just the DBI-action for M5-branes [17] 

Sm5 = -TMS J d^eV-det P[G]a6, (3.9) 
or for M2-branes (see, e.g., [18]) 

Sm2 = -rM2 J d^^V-detP[G]„b. (3.10) 

^■^Our calculations could be extended to fluctuations of the M5-brane in the X'^ direction (parallel to the 
backgroundM2-branes) and of the two-form potential on the M5 worldvolume. Similarly, the M2-brane fluc- 
tuations along the X^- and X'^-directions could also bo considered. In either case, these fluctuations couple 
each other but not to the scalar fluctuations in the 6789-directions. 
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In both cases, P[G\ab denotes the pull-back of the 11-dimensional spacetime metric (3.3) to the 
probe-brane worldvolume. The brane tensions are: tm^ = (27r)~^^p~^ and tm2 = (27r)~^^p~^ 
(see, e.g., [14]). Proceeding as in the ten-dimensional analysis above, we expand these actions 
to quadratic order in the fluctuations to obtain the Lagrangian density 



£ ~ -^/-det gab 



1 + ^^f^'9-^d.X^d,x' 



(3.11) 



where a sum over A is implicit and gab is the induced metric on the probe-brane: 



M5 : ds^ig) = f-^/^{-dt^ + dx\) + /^/^(V + P^d^l + dxj^) , (3.12) 
M2 : ds\g) = -J'^l^df + J^'^{dp^ + p^dO^) . (3.13) 

The induced metric factorizes such that its determinant is independent of the fluctuations. 

Furthermore, retaining terms only to quadratic order in the fluctuations, any dependence on 
the fluctuations can be dropped from the factor f^^^g^'^ in (3.11). Hence once again, the 
quadratic Lagrangian reduces to a simple free scalar theory in a curved background. The 
equation of motion for any one of the fluctuations (^ = 6, 7, 8, 9) then follows from (3.11) 
as 

dc [V-deigabf^^g'^ddx] = • (3.14) 

Now in accord with the above approximations, we take = + and f{r,xn) = 
f{p,xn). It is also convenient to scale the p,xii coordinates by L to define dimensionless 
coordinates g, z: 

p = Lq , xii = Lz . 

Then, apart from an overall factor, the background parameters only appear in the harmonic 
function (B.3) through the ratio L/Ru - see appendix B. 

The uplifted D2-brane solution (r 2> Rn) is applicable when L 2> Rn where we have 

In the field theory, this solution corresponds to the regime gyj^ mq ^ 9ym^^^^ i^]- 
the other hand, the near-core M2-brane solution is relevant for L ^ Rn where 

2^7r^£ 1 

^^'^'^= L6 (l + g2 + ,2)3 - (3-16) 

This solution corresponds to even smaller values of the quark mass in the field theory, i.e., 
ruq <S gyM- absence of the probe branes, the field theory on the worldvolume of 

the M2-branes is super conformal. The probe breaks the conformal invariance introducing an 
energy scale niq. 

We now compute the spectra of mesons corresponding to the transverse scalars for each 
of the probe branes of the uplifted and near-core geometries. Let us begin with the M5-brane 
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probes. For the uplifted D2-brane solution (3.15), the harmonic function depends only on 
and we proceed via separation of variables: 



m 



0,±1,±2,... 



(3.17) 



With = -A;^ and setting 



-2 evr^^^iV , „ 



RuL^ 

the equation of motion (3.14) reduces to the following radial equation for (pup 

M2 £0(^2 + 1) 



(3.18) 



'^up 



dg4>Up{Q) + -dg4>up{Q) 



.(1 + ^2)5/2 

Then the meson mass spectrum is: 

67r2 Nlp^ 



i3 



m 



4>up{q) = 0. 



3 9\m^ 



(3.19) 



(3.20) 



where the dimensionless constants are the eigenvalues of (3.19). In the above, we have 
used eq. (2.9) and the standard formulae (see, e.g.^ [14], [4]): 



(3.21) 



For m = 0, i.e., with no M-theoretic excitations, eq. (3.19) matches precisely the expected 
equation (2.40) for the D2-D4 system and we recover precisely the same spectrum. In the 
M-theory context, one can also excite modes along the xn-circle. Setting $ = (f>up{Q)y^'^ i.S'^) , 
eq. (3.19) can be written as 



(3.22) 



where is the usual Laplacian in three-dimensional spherical polar coordinates. Hence 
eq. (3.22) has the form of a three-dimensional Schroedinger equation with potential V = 
— M^/(l + ^^)^/^ and energy eigenvalue E = —{L/RuYrn?. For small q, the potential 
approaches — while for large it approaches zero as —M'^q"^. Thus, the lowest energy 
eigenvalues will be (very) roughly equal to the minimum of the potential, — M^. Hence we 
would have ~ L?' /R\i and the spectrum of such excitations would be 



m 



(3.23) 



Comparing to eq. (3.20), we note that M^/M^ ~ Lp'/Rii- Recall that the uplifted D2 
geometry applies in the regime L ^ R\\ and so here these m 7^ excitations arc extremely 
heavy and very difficult to excite. As we move towards L ^ Ru, the harmonic function (3.5) 
has a more pronounced dip in the xu direction — see appendix B — and the M-theoretic 
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degrees of freedom become lighter. In the gauge theory, this becomes most pronounced at 
"T-ij 9ym i^-^-t ~ Ru) where instanton effects are unsuppressed [4]. The latter reflect the 
localization of the M2-brane on the xn-circle and so the background goes over to the near 
core solution. 

For the near core geometry (L ^ Rii), the harmonic function for our embedding is given 
in (3.16). In this case, it is useful to use spherical coordinates in the (3,4,5,ll)-space, defining 
the radial coordinate g: 

f = g'^ + z'^. (3.24) 
As usual, we proceed via separation of variables, taking 

X = e"'--y'^{S')cj>core{Q), (3.25) 

where e*'^^ are plane waves in the 01-space and y^^{S^) (^3 = 0,1,2,...) are spherical har- 
monics on the of unit radius satisfying (2.16). With = —k'^ and now setting 

M^= " , (3.26) 

the full equation (3.14) reduces to a radial equation for ^ = (I){q): 

M2 4(4 + 2) 



3 

dl^coreiQ) + -=dg(j)core{Q) + 



(1 + ^')3 ^2 



,iQ)=0. (3.27) 



The mass eigenvalues can again be determined numerically for this equation. The mass 
spectrum of mesons on the codimension-one defect is now 

Clearly, the scaling here has a different form than found in the in eqs. (3.20) and (3.23) in 
the uplifted D2 background. Now, with the identification ruq = L/(27ra') = L Rii/{2Tr£p), 
this mass scale would be ~ mg/{gYM^)- However, we will argue in the discussion 
section that this interpretation is inappropriate and that the correct result, appropriate for 
the superconformal field theory dual to the AdS4 x S""^ core, is ~ ml /N. 

The computation of the spectra in the case of the M2-brane probe is essentially the same 
as the M5-brane case. So in the following we only note the salient differences. In the uplifted 
limit of the background, the ansatz (3.17) is replaced by 

X = e^'^''e^''4>up{6) ■ (3-29) 

The fact that the spatial dependence has been reduced to exp(— zMt) reflects the fact that 
the fundamental fields are localized on a codimension-two defect, i.e., a point-like defect in 
the 2-1-1 dimensions. Scaling the mass as in eq. (3.18), the radial equation (3.19) now becomes 



'-'up 



{g) = 0. (3.30) 
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Hence the spectrum of bound states scales as in eq. (3.20) and in fact, this spectrum exactly 
matches that calculated in section 2.4 with p = 2 since eq. (3.30) is precisely the same as 
eq. (2.63) in this case. Note that since the M2-brane is confined to a point in the xn-circle, 
there are no intrinsically M-theoretic degrees of freedom to be excited at this stage. Rather 
the transverse scalars representing the fluctuations in the .xn-direction simply match the type 
2 gauge modes in eq. (2.70) (while the type 1 modes of eq. (2.71) do not exist for p = 2). Of 
course, this matches the three-dimensional duality relating the gauge field on the D2-brane 
worldvolume to the xu scalar on M2-brane [18]. 

For an M2-brane probing the core limit of the background, the ansatz (3.25) is replaced 

by 

X = e-''^'y'\S^)4>coreiQ). (3.31) 
Scaling the mass as in eq. (3.26), the new radial equation is 



+ 1) 



(1 + ^2)3 



(3.32) 



and one arrives at mass eigenvalues scaling as in eq. (3.28). 



4. Discussion 

In this paper, we used the the gauge/gravity correspondence with flavour to derive the spec- 
trum of mesons for strongly coupled gauge theories in various dimensions. In particular, 
we have studied a Dfc-brane probe inserted into the near horizon geometry of N coincident 
Dp-branes, which in the dual gauge theory corresponds to having introduced a fundamen- 
tal hypermultiplet into the (p + l)-dimensional U{N) super- Yang-Mills theory. The brane 
configuration was arranged to be supersymmetric, still preserving eight of the original six- 
teen supercharges of the background. For each system of Dp- and DA;-branes considered, we 
found that the mesons were deeply bound and the spectra were discrete. Up to numerical 
coefficients, the mass gap for all of these theories had a universal form: 

771 

M ^ , where 5e//("i<?) = 9YMNmP-^ . (4.1) 

Here geff{mq) is the dimensionless coupling constant in the gauge theory evaluated at the 
quark mass niq, which sets the relevant energy scale for the mesons. The full spectrum is given 
as the above mass (4.1) times dimensionless constants, which are computed as eigenvalues of 
an ordinary differential equation. 

These results apply in an intermediate energy regime of the gauge theory where the dual 
ten-dimensional supergravity description is valid. In this regime where g^ff 3> 1, eq. (4.1) 
then dictates that the mesons are deeply bound, i.e., the meson mass is much less than 
(twice) the quark mass. Of course, this result has a simple interpretation from the bulk point 
of view. There, a meson is a bound-state of two strings of opposite orientation (corresponding 
to the quark-antiquark pair) stretching between the probe Dfc-brane and the singularity or 
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horizon at r = 0. To form a meson, the ends at the horizon join together to form an open 
string on the D/c-brane probe. The resulting string is much shorter than the original ones 
corresponding to the 'free quarks', resulting in a lower energy configuration or a meson with 
a mass much less than the quark mass. Since implicitly this discussion depends only on the 
background geometry and the positioning of the probe branc, it is perhaps less surprising 
that the scaling of the meson masses was same for all of the different DA:-branes, even though 
these configurations represent very different gauge theories. 

In fact, we now argue that the universal form (4.1) for the meson masses has a deeper 
connection to holography. In investigating general Dp-brane backgrounds (2.5), one finds that 
probing with wave-packets constructed of massless supergravity fields leads to a 'holographic 
distance-energy' relation [19] which may be expressed as: 

E = , where glff{U) = g^^N U^'^ . (4.2) 

9eff[U) 

Hence the energy of a such wave-packet propagating in the vicinity of (the minimal radius of) 
the Dfc-brane has precisely the same energy as the excitations of the massless open string fields 
on the probe brane. While this may seem a remarkable coincidence, the matching of the open 
and closed string energies is crucial to the consistency of the gauge/gravity duality. Before 
introducing any fundamental fields or probe branes, the duality establishes a connection (4.2) 
between the energy scale in the field theory and the position in the holographic dimension. 
As emphasized in [21], consistency then requires that physics in the two dual theories must 
be local in both of these parameters. In the presence of probe branes, locality of bulk gravity 
theory in the radial direction is, of course, obviously maintained. Locality of gauge theory 
physics in the energy scale was less obvious to begin with, but our conclusions for the meson 
masses (4.1) establish this nontrivial feature is preserved after the introduction of fundamental 
matter. Hence as demanded by holography, we again find the desired locality on both sides 
of the duality. 

We are lead then to conjecture that in fact holography dictates that the energy scale of 
the spectrum of mesons or massless brane fields must be as given in eq. (4.1). This conjecture 
may then be tested in a variety of new contexts. For example, one might consider the 
nonsupersymmetric construction of [7] . There one finds that if the quark mass is much larger 
than the QCD scale (i.e., Ukk) that the meson spectrum takes precisely the desired form. Of 
course, for large niq, the system in [7] is essentially the Dp-D(p -|- 4) configuration of section 
2.2 with p = 4. A more interesting case to consider is the holographic construction of [9] 
which involved a nonsupersymmetric configuration of D4-, D8- and D8-branes. In this case, 
the probe brane configuration is very different from those considered here since the D8-D8 
system connect in a 'wormhole'.^^ However, if one tunes the model, such that the minimal 
radius of the wormhole is above the QCD scale, one again finds that the massive mesons 

^^Rcf. [22] recently considered such probe configurations in a supersymmetric D4-brane background. Ref. [7] 
also discusses a similar wormhole configuration for a D6-D6 pair in a D4 background. 
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have masses scaling as in eq. (4.1) [23, 24].^^ This nonsuper symmetric example, however, 
highlights one small subtlety. The quarks in the dual theory are actually massless, by which 
we mean the 'bare' mass in the UV theory vanishes. However, the infrared dynamics generates 
a finite 'constituent' quark mass, which can be identified with the minimal radius of the D8- 
branes. When we say that the meson spectrum scales as in eq. (4.1) implicitly we are using 
this infrared or constituent quark mass. For a general super symmetric configuration, the 
embedding of the probe brane is flat and so there will be no distinction between the bare and 
constituent quark masses. However, we can expect that in nonsupersymmetric constructions 
that the bare and constituent quark masses will differ as a result of bending of the probe 

— although perhaps not as dramatically as in this example. In such cases, it will be the 
constituent mass that is the relevant energy scale determining the holographic scaling of the 
meson spectrum. As a final comment, we note that the same energy scale appears in phase 
transitions for these systems [25] . 

As indicated above, the full meson spectrum is determined by solving an ordinary dif- 
ferential equation and imposing regularity requirements. The masses are then given by the 
above mass (4.1) times dimensionless constants. While in general one must resort to numerics 
to determine the proportionality constants, it was possible to solve the ODE's (in terms of 
hypergeometric functions) and hence exact values for the proportionality constant for p = 3 

— see appendix A. One of the outstanding features of these p = 3 spectra {i.e., for D3-D3, 
D3-D5, D3-D7) is that they are highly degenerate: all states with the same u = n + £ have 
the same mass. The D3-D7 case was studied extensively in [1] where it was argued that this 
degeneracy may arise because of an extra, hidden 50(5) symmetry. A detailed examination 
of the spectra in section 2 shows that no such degeneracy occurs for p ^ 3 — see, e.g., Table 2. 
Hence it seems that the degeneracy found for p = 3 is connected to the conformal invariance 
of the N=4 super- Yang- Mills gauge theory. 
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21.02 


52.11 


96.90 


155.42 
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8.97 


18.30 


30.51 


45.58 




2 


48.45 


91.32 


147.83 


218.01 






2 


21.47 


35.06 


51.53 


70.89 




3 


86.63 


141.29 


209.53 


291.47 






3 


39.18 


57.01 


77.75 


101.38 



(a) (6) 

Table 2: The spectrum for scalar fluctuations in the (a) D2-D4 and (6) D4-D6 brane conflgurations 
m terms of M^(n,4) . 

The connection between the breaking of conformal invariance and of the degeneracy of 
the meson spectrum can be made more precise as follows. If one considers states with flxed 
u = n + £, one finds in general that the mass increases (decreases) with increasing £ for 

Apart from the massive spectrum, there is also a set of massless pions whose existence is required by 
symmetry-breaking considerations. 
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backgrounds with p = 4 (p < 3) — as can be seen, e.g., in Tables 2a, b. This trend seems to 
be connected to the running of the effective coupUng (1.1) in the nonconformal backgrounds. 
Given that the spectrum scales with M ~ niq/geff, the individual masses will depend on the 
precise value of the effective coupling which is relevant. For a given state, the radial profile 
describes the 'structure' of the meson in energy space and one can calculate the mean energy 
([/) — see the technical details in appendix C. Tables 3a, b show the results for the D2-D4 
and D4-D6 systems, which are of the general case. Note that {U) increases when the radial 
quantum number n is increased but decreases with increasing ^2-^^ Hence combining the 
facts that the effective coupling grows with increasing energy (radius) in the D4 background 
and that the meson masses scale (roughly) inversely with geff, the meson masses would be 
expected to increase more with an increase in i than with a corresponding increase in n. 
Of course, the reverse trend should be expected for the D2 background, where the effective 
coupling decreases at larger energies. Hence this reasoning reproduces the trend in the mass 
spectra commented on above and illustrated in Table 2. This phenomenon is similar to 
the effects seen in the spectroscopy of heavy-quark mesons due to the running of the QCD 
coupling [26]. 
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5.9310 


8.8394 
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2.1634 


3.1276 


4.0868 




2 


1.9985 


2.5708 


3.1427 



(a) 



Table 3: The radial expectation value 2tt {U)/mq for the scalar field x in the (a) D2-D4 and (6) D4-D6 
brane systems. 



Ref. [27] calculated the form factors of the mesons in the D3-D7 system with respect to 
various conserved currents. They found that the size of these bound states to be roughly 
1/M ~ ^ gyjyiN /mq. From the point of view of supergravity, this result is not difficult to 
understand: the form factors are determined by the overlap of various radial profiles and 
M is the only scale in the problem. Hence it is the only dimensionful quantity that can 
set the size of the mesons. This analysis can be extended to mesons in the general Dp-DA; 
systems considered here and the same behaviour will be found with the meson size set by 
1/M ~ gf,ff{mq)/mq. It would be interesting to study these issues in more detail. 

As noted in section 2, we noted that the analysis of the meson spectrum fails for p = 
5 because there are no normalizable solutions. This is related to the usual problematic 
features of holography for D5-branes [4]. At large radius, the coupling becomes large and 



^^The latter is an interesting effect. The 'angular momentum' produces a repulsive, centrifugal barrier, 
causing the peak of the wave function to move away from the origin. However, as £2 is increased, the asymptotic 

potential is also raised causing the asymptotic wavofunction to fall off more quickly. The latter is a stronger 
effect so that the wavefunction localizes at smaller radii and hence smaller (U), as seen in tables 3a, b. 
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the appropriate gravity background is that of the S-dual NS5-brane, where the holographic 
dual is M={l,l) little string theory — see, e.g., [28, 29]. The latter background includes an 
infinite throat containing delta-function normalizable states and so it is natural that analogous 
'meson' states would arise when a probe brane is introduced. This would be an interesting 
topic for further study. 

The computations of the p = 2 spectra in section 3 were motivated by the comparsion with 
previous results for holographic mesons in (2+l)-dimensions [11]. The latter went beyond the 
probe approximation and considered a holographic description of a large number of flavours 
in the strongly coupled D2-D6 system. In particular, their construction began with Nc D2- 
branes and Nf D6-branes alligned as shown in figure (2.30) and used the fully back-reacted 
near-horizon geometry of [20] in the limit 1 <C Nf < Nc but with Nf/Nc finite. In this limit 
with the back-reacted geometry, the fluctuations of the D6-brane correspond to fluctuations in 
the supergravity background and so mesons are described by closed string states in the bulk. 
Ref. [11] attempted to solve for meson masses numerically and found a discrete spectrum with 
a mass gap M oc ruq. The latter seems to contradict the results in section 2.2.2 where we 
found that the meson masses scale as M oc rriq^^ /gy^^- This disagreement is simply resolved 
by realising that the two computations are actually performed for different energy regimes in 
the field theory, as illustrated in figure 1. Our ten-dimensional supergravity result holds in 
the intermediate energy gy^N^^^ <^ niq <^ gy^jN . On the other hand, [11] consider the far 
infrared limit niq <C gyM^ where the gravity description becomes an M-theory configuration, 
i.e., an M2-brane background with orbifold identifications [20]. 

To study meson spectra in different holographic regimes, we turned in section 3 to the 
D2-D4 and D2-D2 systems for which probe brane computations are possible at both the 
intermediate and far infrared energy scales. These correspond to defect theories but the 
results of section 2 show that the mass scale of meson spectra remains the same in these 
cases. In the infrared regime, the D2-brane lifts to an M2-brane and the D4- and D2-brane 
probes lift to M5- and M2-branes, respectively. In the uplifted D2-brane regime, the masses 
of the lowest lying mesons scale as oc L'^Rii/{Ncl^) ~ rriq^'^ /(gyj^^N), which matches 
that found in the ten-dimensional regime. However, we see in eq. (3.28) that this scaling 
changes in the near-core regime: = L^/(iVc^p^). This scaling actually precisely matches 
that found in [11] (once the AdS4 curvature, R\^g = Nclp^, is restored). This agreement is 
almost better than we might have hoped for given that our probe approximation applies for 
Nf/Nc = whereas [11] has Nf/Nc small but finite. 

As explained in the previous section, if we tentatively identified nig = L/(2Tra') = 
L Rn/{2Trip) then the mass scale (3.28) in the near core regime would become ~ nT-q/idyM-^) 
with an extra factor of ruq/gyj^ compared to the intermediate regime where ten-dimensional 
supergravity is valid. However, this is clearly an inappropriate interpretation of the result 
(3.28) which describes the meson spectra in the far infrared where the gauge theory has flowed 
to the superconformal fixed point. In particular, the latter no longer has any knowledge of 
the scales Rn or gvM- The introduction of a probe brane in the AdS4 x 5^ background 
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corresponds to coupling the SCFT to massive fundamental fields, which break the conformal 
invariance (mildly) through their mass. This mass scale must be set by L, the only available 
scale in the supergravity description. Certainly one should not expect that one can extend 
naively extend the identification niq oc L found in ten dimensions to the near core regime. 
In particular, L is now an arbitrary coordinate distance in the core region. Instead wc must 
re-evaluate using the standard distance-energy relation for AdS4 [19], which associates the 
energy with the 'standard' AdS coordinate U. Then with the appropriate coordinate transfor- 
mation, the mass of the fundamental fields in this regime becomes ruq ~ L'^/^p, which matches 
the choice made in [11]. The mass scale in the meson spectrum is then M c± niq/Nc^^. Again 
this result matches the desired 'holographic distance-energy' relation [19]. That is, the scal- 
ing of the meson spectrum here once again conforms to the 'energy locality' demanded by 
holography, as discussed above. 

It is also natural to consider the extension of the p = 4 theories to the strong coupling 
M-theory regime, as well. In this case, the dilaton grows with radius in the D4 background 
and one lifts the the M5-brane background at very high energies. In the dual description, 
this corresponds to the UV completion of the original five-dimensional gauge theory being 
a six-dimensional J\f={0,2) SCFT (on a circle) [30]. If we introduce fundamental fields on a 
codimension two defect as in section 2.4 with p = 4, the D4 probe also lifts to an M5-brane 
(wrapping the eleventh dimension). The fundamental matter excitations are now naturally 
strings associated with M2-branes ending on the probe M5-brane. The low-lying 'mesons' 
again correspond to fluctuations of the massless fields on the M5 probe and one finds that 
the mass gap is again given by iriq/N^^^, in agreement with the 'holographic distance-energy' 
relation [19]. 

As commented above, the agreement of our probe calculations in the AdS4 core of the 
D2-brane background yielded remarkable agreement with the results at finite Nf/Nc in [11]. 
It would be interesting to pursue this further in other backgrounds. Following the results 
of [31], there have been several attempts [32, 33] to establish fully back-reacted backgrounds 
describing four-dimensional gauge theories at finite Nf/Nc- It seems further work is needed to 
develop these solutions to the point where the meson spectra could be calculated. However, 
some simple observations can be made at this stage. In some approaches [33], the back- 
reacted solution still includes flavour branes. However, in other models ek, the branes are 
replaced by deformations of the supergravity background and so the open-string excitations 
representing the mesons would be replaced by closed-string states. On general grounds, 
the latter are expected to satisfy the 'holographic distance-energy' relation in these new 
backgrounds. Hence, from this point of view, it seems that holography will dictate the same 
kind of scaling in the meson spectra at both finite and vanishing Nf/Nc- Hence we expect 
the 'remarkable' agreement observed here for p = 2 might extend to more general situations. 

One direction for future investigation could be considering mesons with higher spin in 
the {p+ l)-dimensional gauge theories. This would require a straightforward extension of the 
calculations for the D3-D7 system in [1]. There the spectrum for large J was calculated from 
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large semiclassical strings hanging down from the probe brane. These spinning string solutions 
can also be extended to consider bound states of quarks with different masses by introducing 
separated probe branes [34] . One universal result will be that in analogy to [1] , the spectrum 
will start with a Regge-like trajectory where the tension is governed by the quark mass. In 
this regime, the gravity description corresponds to spinning strings whose extent is (larger 
than the string scale but) smaller than the curvature scale of the background. The Regge 
slope is then just the redshifted tension of fundamental string and in field theory parameters 
is given by l/ttg^j 2± m^/^e//- While easily understood from the supergravity point of view, 
this behaviour is surprising from a field theory point of view as none of the theories considered 
here are confining [35]. For J ^ 9effi the size of string exceeds the background curvature 
scale. In this regime, it should be that the mass spectrum can be understood as that of 
two nonrelativistic and widely separated quarks which are weakly bound by a long-range 
potential, V ~ —{gyM^/^'^Y^^'^ — the latter being that computed for two static quarks in 
the corresponding background [36] . 
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A. Analytic meson spectra for p — 3 

In section 2 we studied the spectra of mesons corresponding to fluctuations of probe Dfc-branes 
in the near-horizon geometry of N coincident Dp-branes. For each branc configuration and 
each type of fluctuation of the probe brane, we found an ordinary differential equation for the 
radial profile. As discussed in section 2.2, the solutions of these differential equations must 
be real-valued, regular and normalizable in order to be dual to a physical meson state in field 
theory. Thus, solutions were chosen that were real and regular at the origin. The eigenvalues 
M were then determined by requiring that the solutions be convergent as p ^ oo. The 
meson spectrum was then given by (2.19). Though it was not possible to solve the differential 
equations analytically for general p, analytic solutions are possible for p = 3 and we present 
these solutions and the resulting meson spectra in this appendix. 

The D3-D7 brane system was studied in [1] and we review the analytic solutions and 
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meson spectra briefly here. The branes were oriented as shown in the fohowing array: 

0123456789 

L>3 X X X X (A.l) 
Dl xxxxxxxx 

The D7-brane fluctuates in the 89-directions and the radial ODE for these scalar modes is 
given in (2.18). Deflning 

a = -^ + ^v'l + M2>o, (A.2) 

the solution for the radial function (/>(p) is 

4>{p) = p^^ + L^)-" F{-a , -a + 4 + 1 ; 4 + 2 ; -pV^') , (A.3) 

where F{a,b;c;y) is a hypergeometric function satisfying (see, e.g., [37]) 

y(l - y) u"{y) + [c-{a + b+ l)y\ u\y) - abu{y) = 0. (A.4) 

To determine the behaviour of solutions as p ^ oo, the following asymptotic expansion of 
F{a,b;c;y) is useful [37]: 

,. -r^/ , X r(c)r(6 - a) , „ r(c)r(a — 5) , ,1, ^. 

hm F{a,b;c-y) = ( -y « + -^-^ -( (-y) ^ A.5 

T[bjT[c — a) T[ajT{c — b) 

Using this asymptotic expansion, the regularity of (A.3) is achieved by setting —a+i^+l = —n 
(n = 0, 1, 2, ...). The solution is then 

Hp) = (^2 + ^2)n+^3+i n-(.n + 4 + 1) , -n ; 4 + 2 ; -pV^') , (A.6) 
with mass eigenvalues 

=4(n + £3 + l)(n + ^3 + 2). (A.7) 

The solutions for the gauge fields are determined similarly. The ODE for the type 1, 2, and 
3 modes are given in eqs. (2.25), (2.27), and (2.28), respectively. The solutions are 

<^i(p) =p'np' + ^')^''n-« + 4 + l, -a;4 + 2; V/^') , (A.8) 
<^2(p) =p'^-'(p' + i^')~"n-« + 4 + l, -a;4 + 2; -pVl^), 
4>s,+{p) = p'-'+Hp' + L^)-'~''F{-a + 4 + 2 , -1 - a ; ^3 + 2 ; -p^L^) , 
<^3,-(p) = p'^+'(p' + l2)-i-F(-« + 4 , 1 - a ; 4 + 2 ; -pV^') , 



with spectra 



= 4(n + 4 + l)(n + 4 + 2), n>0,4>0 

M| = 4(n + 4 + l)(n + £3 + 2), n>0,4>l 

M|,+ = 4(n + 4 + 2)(n + 4 + 3) , n > , 4 > 1 

M|'_ =4(n + 4)(n + £3 + l), n>0,4>l- 



(A.9) 
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We can see explicitly here that the spectra of the various modes are related as in eq. (2.55). 
As noted in [1], all states with the same n+i^ have the same mass and so the meson spectrum 
has a large degeneracy. 

The D3-D5 brane system was oriented as follows: 

1 23456789 
D3 X X X X (A.IO) 

D5 XXX XXX 

In this system, there are several classes each of scalar and gauge field fluctuations. First, the 

scalar fluctuations of the D5-brane may be parallel to the background D3-branes along the 
direction or orthogonal to the D3-branes in the 789-dircctions. The gauge field modes are 
characterized as type 1, 2 and 3, where the type 3 excitations couple to the parallel scalar 
modes. The ODE for the orthogonal scalars and the type 1 gauge fields are identical and are 
given in eqs. (2.40) and (2.47), respectively. The solutions for these modes are 

0±(p) = 0i(/o) = /'(p' + i^')""^(-a + ^2 + l/2, -a; £2 + 3/2; V/^') , (A.H) 

and the resulting spectra are 

Mi = Mf =4(n + £2 + l/2)(n + £2 + 3/2), n>0, ^2>0. (A.12) 

The type 2 gauge fields satisfy eq. (2.48) and the solutions are 

Mp) = P^'~\P^ + LY"" ^(-« + ^2 + 1/2 , -a ; £2 + 3/2 ; -p^L^) , (A.13) 

with the corresponding spectrum 

M| = 4(n + 4 + l/2)(ra + ^2 + 3/2), n>0, ^2 > 1 • (A.14) 

For £2 = 0, the parallel fluctuations are not coupled to the gauge fields and the ODE reduces 
to eq. (2.49). The solution is 

4>\\{p) = ip^ + LYF{a,a + 5/2; 3/2; -p^/L^), (A.15) 

with masses 

M2 = 4(n + 3/2)(n + 5/2), n > 0. (A.16) 

The type 3 gauge fields and the parallel scalar fluctuations are coupled for £2 > 1- The 
diagonalized differential equations were given in eq. (2.53) with solutions 

4>^(p) = p^^+\p' + l2)-"-i F{-a + £2 + 3/2 , -a - 1 ; £2 + 3/2 ; -p'/L^) 
^_(p) = + L2)-"-iF(-a + ^2 - 1/2 , -a + 1 ; ^2 + 3/2 ; -p'/L^) , 



and spectra 

Ml = 4(n + £2 + 3/2)(n + £2 + 5/2) , n > , ^2 > 1 ; 
Ml = 4(n + £2 - l/2)(n + £2 + 1/2) , n > , £2 > 1 • 



(A.17) 
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The spectra for the various modes satisfy the relations in eq. (2.55) and there is again a large 
degeneracy as the masses only depend on n + ^2- 

Finally, the D3-D3 brane system was oriented as shown in the following array: 

123456789 
background x x x x (A. 18) 

probe XX x x 

Again, the scalar and gauge field fluctuations fall into several classes. The D3-brane probe 
fluctuates in the X"^ and directions parallel to the background branes and also in the 
6789-directions orthogonal to the background branes. The radial ODE's for the orthogonal 
fluctuations (2.63) and for type 1 gauge fields (2.72) were identical. The solutions for these 
fields are 

Mp) = Mp) = + ^')"" ^(-« + 1^1 , -a ; 1^1 + 1 ; -p^L^) , (A.19) 

and with spectra 

Ml = Ml = 4(n + \l\){n + \l\ + l), n > , |^| > 1 . (A.20) 

The diagonalizcd radial equation for ^' corresponding to the coupled X^'^-scalars appears in 
eq. (2.67). The £ > and I < —1 solutions^^ are, respectively, 

<^||,+(p) = /(/ + l2)-"-i F{-a - 1 , -a + £ + 1 ; £ + 1 ; -pVl^), 
<^||,_(p) = p\%p^ + l2)-"-i F{-a + 1 , -a + 1^1 - 1 ; 1^1 + 1 ; -p'/L^) . 

with the corresponding spectra 



=4(n + £+l)(n + £ + 2), n>0, ^ > 0; 
Ml = A(n + \i\-l){n + \i\), n>0,^<-l. 



(A.21) 



Since the ODE defining the spectrum for the ^* modes is eq. (2.67) with I —i, the spectrum 
for these modes is: 



M2_ = 4(n + |^| + l)(n + |£| + 2), n>0, i < 0; 
M^^^ = 4{n + £ - l){n + i) , n>0,i>l. 

For type 2 gauge fields, the radial equation was given in eq. (2.74) and has solutions 



(A.22) 



Mp) = P^'^-\L^ + P')"" ^(-« + 1^1 , -a ; 1^1 + 1 ; -p^L^) , (A.23) 
and the resulting spectrum is 

M'i = A{n + \£\){n + \e\ + l), n>0, \i\ > I . (A.24) 

The spectra of these bosonic modes are related as in eq. (2.75) and as the masses depend only 

on n + 1^1, the spectrum again exhibits a large degeneracy. 

^^Recall from footnote 13 that supersymmetry requires the £ = —1 mode to vanish. 
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B. Localized M2-brane background 

In section 3.1, we studied an M5-brane probe in the near-horizon geometry of N coincident 
M2-brancs. Here, we briefly discuss some details about the near-horizon geometry of the 
M2-branes. 

The metric for N coincident M2-branes is (see, e.g., [14] and references therein) 

where r is a radial coordinate in the 3456789-space and xu is the coordinate in the eleventh 
dimension. We take the eleventh dimension to be compact, with xu = xn + 27ri?ii. For 
simplicity, we consider the harmonic function for the case when all N M2-branes are localized 
at some point along xn — of course, more complicated configurations are possible. The 
harmonic function for the M2-branes can simply be written as a sum over images in the 
xii-direction [5]: 

It is relatively straightforward to sum this expression in closed form using techniques from 
complex analysis (see, e.g., [38]). The final result can be written as: 



f(,r, Xll) 



p_ 



Rii sinh(r/iiii) o-^ii cosh(r/i?ii) cos(xii/i?ii) — 1 



+ 3 



cosh(r/i?ii) — cos(a;ii/i?ii) (cosh(r/i?ii) — cos(xii/i?ii))^ 

sin(r/i?ii)(2 - cos^(a;ii/i?ii) - cosh(r/i?ii) cos(a;ii/i?ii))l 

(B.3) 



(cosh(r/i?ii) — cos(xii/i?ii))3 



As in the fluctuation analysis in section 3, we take = + and scale the p and xn 
coordinates by L: 

p = Lg , xu = Lz , 

where the z-periodicity is given by z ~ z + 27ri?ii/L. The harmonic function (B.3) can then 
be expressed as 



27r^^/iV 



],R'u 1 sinh(^yrT7) 

L5 (1 + ^2)5/2 ^^gj^(_L_yi^))_^^g(_L_^) 



+3^1 1 cosh(^7rT7))cos(^z)-l 
(1 + (cosh(^yrT^)) - cos(7^.))2 
Rl^ 1 sin(^yr + 7))(2 - cos^(^^) - cosh(^7rT7)) cos(^z))^ 

(1 + ^2)3/2 (cosh(^yrT^)) - C0S(^Z))3 

In section 3.1 we focused on the regimes L Ru and L <^ Rii, which we referred to 
as the "uplifted D2-brane" and "near core" solutions, respectively. Figure 2 displays plots 
of the full function (B.4) for various values of L/Rn to illustrate how the full harmonic 
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function changes in different parameter regimes. Figure 2a displays the harmonic function 
for L/Rii = 50. This ihustrates how in the L ^ Rn regime the structure in the xn 
direction is washed out and / essentially depends only on the radial coordinate p. That is, 
the full harmonic function is well-approximated by the uplifted D2-brane solution (3.15). The 
harmonic function in the limit L <ti Rn is shown in figure 2c, where L/Ru = 10"^. In this 
case, the structure in the p and z coordinates is the same and f{p,z) can be approximated 
by the near core solution (3.7). Figure 2b shows the harmonic function for L/Ru = 3, an 
intermediate regime. 



Figure 2: The harmonic function —Rii/(2TT^£p^N)f{p,z) (where f{p,z) is given by (B.4)) for (a) 
L/Rii = 50, (b) L/Rii = 3, and (c) L/Rn = 10""^. 

C. Radial expectation values 

In order to compute the expectation value (r) (of the radial coordinate), the radial eigenfunc- 
tions must be normalized. Here, we normalize the modes following the prescription of ref. [39]. 
We briefly review their formalism (which is a simple generalization of the usual normalization 
of modes in four dimensional quantum field theory) before modifying it appropriately for our 
uses. 

Consider a scalar field propagating in a (c? + 1 + A;) -dimensional spacetime with metric 




rho 



(a) 



(b) 



(c) 



ds^ig) = f{y)v,iudx^'dx''+g±rnndy"'dy 



,71 



where p^u = ...d, m,n = 1, k. For a scalar field with the Lagrangian density 





the canonical commutator is 



$(x,y),6(x',y') 
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Specializing for the moment to the case d = 4 and expanding the field as 

with the oscillators covariantly normalized as [aa(k),aj(k')] = {2TrfS(^\k-]<.')Sa/3, the modes 
in y-space are then normalized according to 



/ 



d'y^MvW^iy) = Sap- 



f 

For any one of the Dp-Dk brane configurations, the Lagrangian density for the scalar fields 
X corresponding to fluctuations of the D/c-brane transverse to the Dp-branes was (omitting 
multiplicative constants) 

£ = e-^^~dctga,J-'g%daXdbX (C.2) 

where gab is the induced metric on the DA;-brane. Comparing (C.l) and (C.2) suggests that 
for the scalar field confined to the D/c-brane, the factor y/— det gab should be replaced by 
e~'^\/— det gab f"^- Thus, the canonical commutator should be 

lx(^,yU«v')] = «pJ=pr*^. - .')sHy - V') 

and with x expanded in terms of Fourier modes and spherical harmonics (as described in the 
text), the measure for the normalization for the radial functions (t){Q) should be 

I^{q) = e-V-det^ab/-'. (C.3) 
Hence, the radial expectation value of any function will be given by 

where the radial measure for any one of the Dp-Dfe brane systems is 



.1+ 



k — p 



Given the standard duality between energy and radius, i.e., U = r/a', the radial pro- 
files give the wave- function in energy — that is, ^(£))| </>(£)) P serves somewhat like a parton 
distribution function. The expectation value (r) gives a measure of the mean energy of the 
components of a given meson state. Using (C.4) and (j){Q) for the transverse scalar fluctu- 
ations of the probe brane, we computed 2n{U)/mq = {(1 -|- Q^Y^'^) for the T>p-T)k systems. 
Note that for p = 3, the calculations can be done analytically while for p 7^ 3, the results are 
obtained numerically. Explicit results for the D2-D4 and D4-D6 configurations are displayed 
in Table 3. These are typical of the general case. Note that {U) increases when the radial 
quantum number is increased, but decreases with increasing £2- 
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In [1], it was noted for the D3-D7 system that for large R-charge, £ ^ n, the radial profiles 
were sharply peaked at Qpeak = 1- One can see that the same peaking arises for the mesons on 
other probe branes in the D3 background from the analytic solutions, e.g., , eqs. (A. 8), (A. 11) 
and (A. 19). This radius corresponds to the radius of a mill geodesic on the induced probe 
brane metric orbiting on the internal space at a constant value of g. One expects a similar 
peaking for occurs for general for wliicli the geodesic analysis yields Q^qq^. — 

2/(5 -p). For 

the large R-charge limit, £ 3> n, one can then infer 2'k(JJ) /vriq ~ {7 — p) / {b — p). 
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